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Abstract 

We give a self-contained presentation of a novel approach to a construction of spherical 
harmonic expansions on the unit sphere in C". We derive a new formula for coefficients of the 
expansion of a smooth zonal function defined on the unit sphere and apply it in some special 
cases. The expansion for the Poisson-Szego kernel for the unit ball in C n obtained by our 
method coincides with the result obtained originally by G. Folland, and on the other hand 
disproves results recently presented in a paper of V.A. Menegatto et al.. 

AMS Math. Subj. Class. (2000): 58G35, 35F05, 33A75, 33A45 

Key words and phrases: Laplace operator; spherical harmonics; zonal harmonic polynomials; 
Jacobi polynomials; Poisson-Szego kernel. 

1 Introduction 

This paper presents a part of results from the doctoral dissertation of Agata Bezubik [2] concerned 
with spherical expansions of zonal functions on the unit sphere in the complex space C n . The aim 
was to develop further an approach initiated in the papers [3, 4] of constructing spherical expansions 
of smooth zonal functions in terms of their differential characteristic (e.g. the Taylor coefficients) 
rather than the integral ones (i.e. integral means against appropriate reproducing kernels). 

The complex case treated here, in contrast to the real case discussed in the earlier paper [4], 
goes beyond the classical theory in the sense that it requires the use of an orthogonal system in two 
variables — more specifically, of an orthogonal system on the unit disc constructed in terms of the 
so called disc polynomials obtained from the classical Jacobi polynomials, cf. e.g. [10, 15]. However 
this aspect is interesting for us only so far as being the tool for parametrizing the zonal functions 
on the unit sphere in C n . 

Finally let us note that almost at the same time as the doctoral thesis of A. B. was submitted 
(June 2010), an article [11] of V.A. Menegatto, A. P. Peron, and CP. Oliveira has been published 
in this Journal (Collectanea Mathematica) . Acknowledging an inspiration of our earlier paper [4], 
the authors formulate an expansion theorem (Theorem 2.3 of their paper) with an almost identical 
content as the Theorem 2 — this ,, almost" referring to certain inaccuracies regarding the various 
coefficients throughout the paper. We shall pinpoint these problematic expressions in [11] in due 
course below. 
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2 Preliminaries 

We shall be concerned with problems in analysis on the unit sphere S = {z £ C | |z| 2 = 1} in the 
complex n-space C™. The norm | • | is the usual one, defined as \z\ 2 — Y^j=i \ z j\ 2 ' an d derived from 
the hermitian inner product (z\w) — Y^j=i z j^j- Given z = (z\, . . . , z n ) £ C™, we set zj = xj +iyj, 
with real Xj, yj and i = for j = 1, . . . , n. This gives us a standard identification of C™ with 

the real Cartesian space R 2 ™ by (zi, . . . , z n ) -n- {x\, . . . , x n , yi, . . . , y n ). 

The special unitary group SU(n) acts transitively on S with the isotropy group of a point tjeS 
isomorphic to the group SU(n — 1). In fact, taking r\ — e„, the unit vector corresponding to the 
ra-th coordinate, we easily see that the isotropy group consists of matrices of the form ( q \ \ with 
the (n — l)-square block A in the upper left corner belonging to SU(n — 1). 

Thus, along with the usual representation of S as the homogeneous space S — SO(2n) / SO(2n — 1) 
of the orthogonal group SO(2n), we can identify the sphere S with the homogeneous space SU(n) / SU(n — 1). 
The choice of reference point i] e S is not relevant, since isotropy groups are conjugate and it will 
be of some advantage not to fix this point by any specific choice. The sphere S is equiped with the 
normalized Euclidean surface measure do2n-\ = da (i.e. J s da = 1), which is 0(2n)-invariant, as 
well as SU(n)-invariant. 

By means of the above identification we can express objects attached to the Euclidean space R 2n 
in terms of the complex coordinates. For example the Euclidean Laplacian A on R 2 ™ is written in 
the complex form 

n f) 2 

using the standard notation B- = \(jL.- , = \ (jL. + i^j for j = 1, . . . , n. 

Likewise, polynomials on R 2 " can be written as polynomials of the complex coordinates Zj and 
their conjugates z~j, what will be indicated by writting P(x,y) = P(z, z) for a typical element of 
the polynomial algebra ^(R 2 ™) also denoted as "P(C"). For example, the square of the norm is 

A departure point for analysis on the unit sphere S d 1 is the fundamental decomposition 

oo 

L 2 (S , da) = 0# z . (2) 

1=0 

where H l denotes the space of restrictions to the sphere of harmonic (i.e. annihilated by the 
Laplacian) and homogeneous of degree I polynomials on R d . A key property of those spaces is their 
irreducibility under the action of the orthogonal group SO(d). 

For the case of the unit sphere in C™, the main role is assigned to the group SU(n), what requires 
a refinement of the decomposition (2) based on the notion of bihomogeneity of polynomials. It is 
introduced as follows. Given a pair of non-negative integers p, q one says that a polynomial P(z } z) 
on C" is bihomogeneous of degree (p, q) if 

P(Xz, Xz) = \ p X q P(z,z), AeC, zeC". (3) 

The space of bihomogeneous of degree (p, q) polynomials on C" is denoted by V M = T {p ' q) {C n ). 
We notice that P^p^ contains only holomorphic polynomials, while p(°' q } — antiholomorphic ones. 
Clearly, (p, q) bihomogeneous polynomials are homogeneous of degree p + q in the real sense and 
there is a direct sum decomposition 

P l (C n ) = © -p(P>9)((C rl ) 
p+q=l 

of the space of homogeneous of degree I polynomials on C". We shall write = H( p > q \C n ) = 

ker AnV ( - p ' q \C n ) for the space of harmonic and bihomogeneous polynomials of bidegree (p, q) (with 
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C™ as their domain) — termed " solid harmonics" , and H^ p,q ^ for the space of their restrictions to S 
- "surface harmonics". These spaces are irreducible under the action of SU(n), their dimensions 
are given by (c.f. [10, 15]) 

dimW<™> = dimtf<™> = (n+p + q -l)( n -2+p)l( n -2 + q y 

p\ q\ (n- 1)! (n- 2)! 

and there is an orthogonal decomposition (w.r. to the L 2 (S) inner product) 

H l = (4) 

p+q=l 

and consequently a SU(n)-irrcducible, orthogonal decomposition 

oo 

L 2 (S, da) = © (5) 

Analogous to the real case, there is a decomposition of bihomogeneous polynomials into harmonic 
components. In fact, A : (C n ) 7?(p-i>9-i)(C") is surjcctive with the kernel "H (p < 9) (C") and 

consequently there is a SU(n)-invariant decomposition 

min(p,g) 

V {p ' q) (C n )= r 2fc -H (p ~ M ~ fe) (C"). (6) 

fc=0 

The explicit formulae for this decomposition can be deduced from their counterparts in the real case, 
which were stated in our paper [4, Theorem 1]. We just formulate the end result below. It should 
be noted, however, that the decomposition has a long history, as it was studied independently by 
several authors, see e.g. [9, 17, 10]. 

Theorem 1 (The canonical decomposition of bihomogeneous polynomials) Given p,q e 
Z + we set m = min(p, q) and for < k < m we define 



p, q ), ~ = . 1 v (n - 1 + P + q - 2fc)(n - 2 + p + q - 2k - j)\ 
k (J) [ ' A k +3 k\j\ {n-l+p + q-k)\ 



and for a polynomial P e "P*^ 1 ?) we se t 

m — k 

h ( -r\p) = j2^' q) (jy j ^ +j (p)- (7) 

Then h { £' q) {P) G %(p- k 'i- k ) and P can be decomposed into harmonic components 

m 

p=j2r 2k h { r\n (8) 

The resulting direct sum decomposition ptv^ = Q)™ =0 r 2k 'H.( p ~ k ' q ~ k ) is orthogonal with respect to 
the inner product induced by restriction of the inner product in L 2 (S) and the maps 
r 2k h^' q \p) are orthogonal projections corresponding to this decomposition and commuting with the 
action of SU(n). 
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3 Decomposition of zonal functions 
3.1 Preliminaries 

Recall that the isotropy group SU(n) of any given point r/ E S = S 2n ^ 1 C C™ is conjugate to 
SU(n- 1). A function / defined on S which is invariant under the action of the isotropy group 
SU(n) will be called a zonal function with the pole at rj (more properly perhaps it should be 
called "zonal with respect to SU(n)"). By a slight abuse of notation we shall call a bihomogeneous 
polynomial on C" zonal, if its restriction to the unit sphere is zonal. Since the space SU(n) \S of 
orbits of the action of SU(n) on the sphere S is naturally identified with the closed disc D — {w € 
C : \w\ < 1} in the complex plane by factoring the map 

S 3 ^ (£ I rj) e D (9) 

with respect to the natural projection S — > SU(n) ?) \S, it follows that any zonal function / on S 
can be uniquely represented in the form 

m = m\v)), t^s. (io) 

by a function (p defined on the disc D, called the profile of /. 

Further, the (normalized) measure dazn-i on the sphere S is factorized with respect to the map 
(9) as follows 

/ /(Od*2n-l(0 = — [ (l-\w\ 2 ) n - 2 d\(w) [ f(wT,+ y/l-\w\*p)d<T 2n _ 3 (p), (11) 

where d\(w) is the Lebesgue measure on the disc D and S 2n ~ 3 is the unit sphere in the space 
perpendicular to rj. 

To describe zonal functions we need to introduce a convenient notation and recall some classical 
notions. We write a = n — 2, since in this form the dimension of C" enters most of the formulas 
in the sequel. The hypergeometric function is defined by the sum of the series, convergent for all 
1*1 <1, 

2Fl U ' V-^-wTi!' (12) 

where a, b, c are parameters (for us always real), and (r)j denotes the Pochhammer symbol defined 
recursively by setting (r)o = 1 and (r)j = r(r + 1) • • • (r+j — 1), for j > 0, for any real (or complex) r 
and non- negative integer j. In the following we shall be concerned with the case when the nominator 
parametres a, b are non-positive integers and c < min(a, 6), in which case the hypergeometric series 
(12) terminates and represents a polynomial of degree min(— a, —6). 

The next result provides at the same time explicit expressions for zonal harmonics and the 
detailed form of the harmonic decomposition of certain bihomogeneous polynomials. 

Proposition 1 Let p, q £ Z + be given and a = n — 2 as defined above, fix rj £ S and consider the 
zonal polynomial P(z,z) = (z \ rj) p (z \ rj) q € "P^'^. The harmonic decomposition (8) of P has the 
form 



(z\rjr(z\rj) q =J2\z\ 2K 



p\ q\ (a + l+p + q-2k)\ 



k\{p-k)\(q-k)\(a + l+p + q-k)\ 
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The proof is a straightforward computation based on the formula (7) combined with some chasing 
of coefficients and can be omitted. 

To note is the representation theoretic interpretation of the decomposition (13). Since the 
projection maps P i— > h^' q \p) commute with the action of SU(n), the harmonic polynomials 

/ , \p-fc7 — i — Ti- k t? ( -P + k, —q + k 1 \ 

in the above decomposition are zonal elements of SU(n)-irreducible spaces %(p- k ^^ k ) . Since the 
latter are known to contain one dimensional subspace of zonal elements, we may state the following. 

Corollary 1.1 Given arbitrary nonnegative integers k, j, any zonal polynomial in %( fe >-?) with pole 
at rj is proportional to the following 

4 fe '%- *) - (* i ntvwi^ (_- k : k ] - r k^rf) • (14) 

It is customary to normalize the zonal polynomials by requiring that they assume the value 1 at the 
pole. Since the Chu-Vandermonde formula (cf. [1, Cor. 2.2.3, p. 67]) gives 

iFi ( ~k, ~j . x \ _ {a + l)k{a + l) 3 



a-k-f J {a + l) k+j ' 

we set 

obtaining so called reproducing kernels for spaces H^ k '^. They are characterized by the identity 

dim / f(OZi k,j) (OMO = f(v), f G ^ (fe,j) - (16) 
Js 

An equivalent expression for zonal kernels uses Jacobi polynomials Pm' v \t), cf. [5, Section 2.4.3], 

-Lm (t j 

where m = min(fc, j). Taking account of the normalization used, we can rewrite the decomposition 
(13) for surface harmonics as follows 



7? _ 

k=0 

with 



_ p\ql(a + l+p + q-2k){a+p-k)l(a + q-k)\ 
lk k\a\(a + 1 + p + q - k)\{p - k)\{q - k)\ ' [ ' 

Profile functions of the reproducing kernels Z^'^ (£) are precisely the disc polynomials in the 
normalization given in the book of Dunkl and Xu [5, Section 2.4.3]. For p, q e Z + and a as above 
we set 

= 7 ^it!^^ wP ^ 9 ^ ( ~ P ' ~ q : A V w e A (19) 
p ' 9 (a + l) p (a + l)g ^-a-p-g' |w| 2 / 

so thatZ^« ) (0 = W p %((» 7 |0)- 

In virtue of surjectivity of the map (9) the decomposition formula (13) implies the following 
decomposition of basic monomials on the disc. 
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Corollary 1.2 For any pair of nonnegative integers p, q, the following is valid 

m 

Ww" = J2^' 9 W^ kiq _ k (w), weD, (20) 

fc=0 

where 7^' 9 are given by (18). 

Substitution of w = 1 in (20) together with the normalization W^ p ^ (w) — 1 gives the following 
"decomposition of unity" 

_ ^ p\q\(n -l+p + q- 2k) {n -2+p- k)\(n - 2 + q - k)\ 

~ £^ k\(n - 2)\(n - 1 + p + q - k)\(p - k)\{q - k)\ ' ( ' 

which is crucial for the proof of our main result, the Expansion Theorem 2 in the next Section. 

Remark 1 The expansion (20) coincides, with the necessary adjustment of notation, with the one 
given in the paper of Wunsche [18, formula (3.10), p. 142], where it is approached from a different 
angle. However, in the paper of Menegatto and al, [11, Lemma 2.1, p. 153] Wiinsche's formula 
is incorrectly quoted, with an extra factorial appearing in the expression for their coefficients 
(replacing our^ q ). As a consequence, their decomposition of unity (2.5) on page 154 does not 
agree with our formula (21), and is not correct, as the following simple example shows. Taking 
m = q = 2,n = 3 in their formula (2.5), one gets 

mi y^' 3) 2!3!(2 - k + 2 - 2)!(3 - fc + 2 - 2)!(2 + 3 - 2k + 2 - 1)! _ 
f^ Q (2 - 2)!fc!(2 - fc)!(3 - fc)!(2 + 3 - k + 2 - 1)! ~ ' ^ ' 



3.2 A general expansion theorem 

We now come to the main result of this paper, which generalizes Theorem 2 of our earlier paper [4] 
to the present context of spheres in complex n-space. 

Let / be a zonal function on the unit sphere S C C™. Recalling this means that / is invariant 
with respect to the isotropy group SU(n) of a point rj e S, it is easy to infer from the group 
invariance of the decomposition (5), that its spherical harmonic expansion is of the form 

00 

f(0= E d P^ dim 4^)(£), 

p,q=0 

where d Pj q — d Ptq (f) are scalar coefficients. Now, regarding the issue of computing these coefficients, 
the general theory of orthogonal expansions together with the identity (16) implies that 

d P , g (f)= I fm^ q) (0da{0- (22) 
Js 

By virtue of (11) this reduces to the integral of the profile ip of /, 

d P , q (f) = — / <f(w)WZ q (w)(l - \w\ 2 ) a d\(w). 

Similarly to the case of the real sphere, discussed in our earlier paper [4], for sufficiently regular 
zonal functions the integral formula for the expansion coefficients can be replaced by a differential 
formula involving Taylor coefficients of the profile function. 

To simplify notation we shall write d, 8 instead of , ^ . 
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Theorem 2 Let f be a zonal function with pole at rj £ S, <p, /(£) = | if)). Assume the profile 
function <p is real analytic on the interior of the disc D and its Taylor series around 



E 



y(0) • k 



j,k=0 J 

is absolutely convergent on the unit circle. Then the coefficients d Pt q of the spherical harmonic 
expansion of f, 

oo 

/(O = E d P>i dim^fe") Z^)®, (23) 

p,q=0 

are given by the formulae 

^-^'-"- l » i |tH» < ^CC»)r (24) 

and the expansion is absolutely and uniformly convergent on S. 

Proof. The assumption on the Taylor series of (p implies that the series converges uniformly to ip 
on the closed disc D, hence we can write 

~ a^V(o) 7 -_ fc 



j\k\ 

j,k=0 J 



Using the formula (20) we get 



j,k=0 J ' 1=0 



Substituting w = (£ | rf) turns the above into 

/«) - E ^^E'T/ , *^- ,lfc - ,, (0, - - nhn(,, k). (25) 
i,fc=o ■ l=0 

To justify the change of the summation order it is enough to show that the inner sum is uniformly 
bounded (w.r. to j, k). Since 

| Z ti-i,fc-D ( £)| = |W7-;, Ml <1, 
cf. [5, 2.4.3, Property (m), p. 57], we get 



z=o 



<E7r fc = i. 



using (21) at the end. 

Now in the expansion (25) we first sum the terms proportional to the zonal harmonic zjf' q \^) 
of a given bidegree. The resulting coefficient is 

- / aH-'fr+y ( o) + , g+ A (P;9) 
2-.{( p + iy.( q + iy l1 ) n w 
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Using (18) we get 

1 p+Lg+i = (a + l+p + q)(a + p)!(a + g)! 

{p + l)\(q + l)\ l1 p\qWa\(a + l+p + q + iy. 

= dim H<*^ + 

l\(a + l+p + q + l)\ 

This finishes the proof. □ 

Remark 2 This result was obtained by one of the authors [A.B.] of the present paper in the course of 
the work on her doctoral dissertation, which was submitted to the Warsaw University of Technology 
in June 2010. An analogous result is given in the paper [11] of Menegatto and al., however the proof 
given there rests on the use of their decomposition of unity (2.5) on page 154, which is not correct 
as it stands. 

3.3 A variant of the Funck— Hecke formula 

The expansion formulae (23-22) immediately imply the following form of the classical Funk-Hecke 
theorem for complex spheres, cf. [13, Theo. 4.4]. 

Corollary 1.3 If f is a zonal function /(£) = </?((£ |?y)) satisfying the assumptions of the Theorem 
2 above, then for every complex spherical harmonic Y( p > q ) e H^ p ' q \ 

f Y^"\Of(0 MO = d p , q ^)Y^ q \ v ), (26) 
Js 

where 

Proof. In fact, substituting the expansion formula (23) into the integral on the left-hand-side of 
equation (26) and changing the order of integration and summation we get 

/ Y^(0f(0do-(0= £ d p , q dim f Y^"\0^ l9 H0M0 

By virtue of (16) and the orthogonality of spherical harmonic of different bihomogeneity we get the 
result. □ 

Comparing with the formula (4.3) of Quinto, [13, p. 258], see also (11) above, we see that 

d p+k d q+k ip(0) 



W^> q Xw)y{w){l - H 2 )™- 2 d\(w) = n(n - 2)!^ - 



This is reminiscent of the classical Pizzetti's formula, cf. [12], however this connection will be 
considered elsewhere. 



3.4 Applications 

In the following we give two applications of the main theorem. 
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3.4.1 Plane wave expansion 

Corollary 1.4 Given rj € S let x i-> e lRc ( x l')) he the plane wave in C n with normal r\. Then writing 
x = \x\£ with £ G S, one has the complex spherical expansion 

e iRe ^) = (n-l)!(^i) " ^ 9 <^HMj^ q+n -i(\x\) ZM(S), (28) 

p,q=0 



where the Bessel functions J u {r) of the first kind and order v are given by 

r(k + - 

The expansion (28) is absolutely and uniformly convergent on every ball in C r ' 



*<"-(£) giWWTTTTT)© ' <29) 



Remark 3 The Bessel functions, which are coefficients in this expansion, depend only on the total 
degree of homogeneity p + q. Therefore the expansion (28) can be reduced to the usual plane wave 
expansion in R 2n , cf. eg. [7, 3} 

e iR ^ = (n-l)\(Jf) " ^dim^+^M)^) 



z=o 

by using the summation formula for the zonal kernels 



dimH 1 4(0 - Yl dim ^ (p,9) z n"' q) (0- 



p+q=l 



Proof. By differentiation of the profile function (f(w) = e lrKcw , where r = \x\, we obtain for the 
coefficients of the plane wave expansion the following expressions 

°L (_Y\k -p+q p+q+2k 

d p , q (<p) = (n-iy.^ 



2P +( i+ 2k k\(n -l+p + q + k)\ 

(n i)< v+«(LY n+1 (LY +9+n ~ 1 V til! (L) 

~ [H \2) \2) k\{n-l+p + q + k)\\2) 



fc=0 

By comparison with the formula (29) the result follows. □ 
3.4.2 Expansion of Poisson Szego kernel 

Let us remind that the Poisson-Szego kernel is defined by the formula 

fl - \z\ 2 \ n 

where B C C™ is the unit ball and the unit sphere S is its boundary. It plays the same role with 
respect to the Laplace — Beltrami operator associated to the Bergman metric on B (cf. eg. [14]) as 
the usual Poisson kernel in relation to the Euclidean Laplacian in the real case, i.e. the formula 



u(z) = / P„(z, 
Js 



V)f(v)da, zeB, 



expresses functions in the unit ball annihilated by the Laplace-Beltrami operator in terms of their 
boundary values on the sphere S. 
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In the paper [8] Folland derived the spherical harmonic expansion of the Poisson-Szego kernel 
and gave the explicit expressions for the coefficients in terms of hypergeometric functions. We show 
below that Folland's expansion follows from our expansion Theorem 2, which applies here in view 

(1 - r 2 )" 

of the fact that the function £ i-> P n (r£, ri) = -. j-. is a zonal function on S with the pole 

|i-K£h)l n 

at n. 

Theorem 3 For r e [0, 1) and r/ e S the spherical harmonic expansion of the Poisson-Szego kernel 
£ i y P n {r£, n) has the form 

P nH,v) = u_ r{ lL n = E dim^SrW^CO, (31) 

I VS I // I p g—Q 

where 

cp.q (r) = rP+q (p + n-m q + n-l)\ ~ (p + k - l)\(q + k - 1)\ r^ 
n{> (n-l)!(p- 1)!( 9 -1)! ^ (n-l+ P + q + k)\ k\ 1 ' 

or equivalently 



n[ > \p + q + n' J J 21 \p + q + n' J 



(33) 



Remark 4 Folland in [8] approaches the problem of this expansion along different lines. He first 
establishes a special case of the solution of the Dirichlet problem corresponding to a bihomogeneous 
spherical harmonic as the boundary value — at this stage already the series S^' 9 (r) for the coefficients 
enter, and from that derives the expansion of the kernel by examining the convergence of the right 
hand side of (31). 

Proof. From the expansion given as formula (2.9) on p. 6 in [10], cf. also [15] 



E 



p + n — 1\ f q + n — 1 



\l-w\ 2n (l-w) n (l-W) n p Z ^ Q \ n-1 J V n ~ l 
setting ip r {w) = |1 — rw\~ 2n we immediately see that 

d>d k <p r (P) =r ]+k {n) 3 {n) k . 

Now applying formula (24) of the Theorem 2 to the function t/?r((Ch)) (with rj <G S kept fixed) we 
see that the following holds 

Lemma 2 For any £, n e S 



l-r(£\v)\ 2n 

- 

Jn~- l)!(n- 1+p + q + ky. k\ 



V dim Z<™> (£)r^ V (n- 1 + p + k)\(n - 1 + g + k)\ r™ 



p,9=0 fc=0 

To finish the proof of Theorem 3 we need to expand( 1 ) 

(n-1 +p + k)l(n -l + q + k)\ r 2k 
; (i (n-l)\(n-l+p + q + k)\ k\ 



(1 -r 2 )™r p+9 E 



in terms of powers of r 2 and show that the expansion coincides with the formula (32). 



1 There is no attempt to make this expansion in the paper of Menegatto ct al.[ll] 
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Using the binomial formula we get 

(n - 1 + p + k)\(n - 1 + q + k)\ r 2k 
(n-l)\(n-l+p + q + k)\ k\ 



oo 



(1 -r 2 ) n r p+q J2' 



k=0 



n / \ oo 

£ ( -iW"W 

'■ -i— n J.— n 



(n — 1 +p + £;)!(« — l + (j + fc)!r 



I Jlk 



{n-l)\j^ \jj {n-l+p + q + k)\ k\ 

rP+Q V r 2 fc y, v fn\ (n-l+p + k-jy.(n-l + q + k-jy. 

^ Hi/ (n-l+p + g + A-i)!(fc-i)! 1 J 

The inner sum in the last equality, after rewriting it in the form of a hypergeometric sum, can 
be computed by using the Pfaff-Saalschutz identity, cf. [1, Theorem 2.2.6], 

( -n, a, b \_ " (-n)j(a)j(b)j 1 _ (c-a) n (c-b) n 

3 2 \c, 1 + a + b-n-c ' J (c)j(l + a + b - c - n)j j\ (c)„(c - a - b) n ' 



In fact, recalling the identity satisfied by the Pochhammer symbol (— l)- 5 — — ^7 = (— t)j, for 



t > j, we see that the inner sum in the expression (35) can be rewritten as a hypergeometric type 
sum 

(n- l+p + k)\(n- 1 + g + fc)! (-r^ji-tyjjl-p-q-k-n), 1 

k\{n- l+p + q + k)\ (l-p-k-n)j(l-q-k-n)j j\' [ ' 

Thus setting a — —k, b=l— p — q — k — n, c=l— p — k — n'm the Pfaff-Saalschutz identity and 
using (t) n = (— 1)™(1 — t — n) n , the sum (36) is seen to to be equal to 

(p + n- \)\{q + n- l)\(p + k - \)\{q + k - 1)! 
k\(p - l)\(q - l)\(p + q + n + k - 1)\ 

Substituting this expression into (35) and comparing with (32) concludes the proof of the Theorem 
3. □ 

Remark 5 The expansion of the Poisson-Szego kernel given in the Corollary 3.3 of [11] differs 
at two important points from the expansion (32), which coincides with the one given originally by 
Folland in [8J. Firstly, on the right hand side of the given formula the factor (1 — r 2 ) q is not 
expanded, but more importantly, the inner sum there is just a numerical factor, because the factor 
r 2j is missing in all its terms. The consequence of this omission is that the formula (3.11) in the 
Theorem 3.4 of the paper [11] expresses S^' g (r) as a polynomial, what is definitely not true. 
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